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The effect of proximity to a Mott insulating phase on the charge transport properties of a super- 
conductor is determined. An action describing the low energy physics is formulated and different 
scenarios for the approach to the Mott phase are distinguished by different variation with doping of 
the parameters in the action. A crucial issue is found to be the doping dependence of the quasipar- 
ticle charge which is defined here and which controls the temperature and field dependence of the 
electromagnetic response functions. Presently available data on high-T c superconductors are anal- 
ysed. The data, while neither complete nor entirely consistent, suggest that neither the quasiparticle 
velocity nor the quasiparticle charge vanish as the Mott phase is approached, in contradiction to 
the predictions of several widely studied theories of lightly doped Mott insulators. Implications of 
the results for the structure of vortices in high-T c superconductors are determined. 



I. INTRODUCTION 



High-r c superconductors are created by doping an antiferromagnetic 'Mott insulating' parent materialEl. A Mott 
insulator is a material which band theory would predict to be metallic but which, because of electron-electron interac- 
tions, is insulatingp. It is conventional to distinguish the different doping levels as 'overdoped', 'optimally doped', and 
'underdoped'. Optimal doping is generally defined as the carrier concentration which maximizes the resistively-defined 
superconducting transition temperature T c . Overdoped materials have more added carriers and underdoped materials 
have fewer. All high-T c materials exhibit behavior which deviates from the 'Fermi-liquid-Migdal-Eliashberg' theory 
which describes most metals; however, the deviations grow more pronounced as the doping is reduced towards the 
Mott insulator, and indeed understanding the physics of the underdoped regime (in other words, of the approach to 
the Mott insulating phase) has emerged as one of the crucial questions in the high-T c field. 

High-r c superconductivity remains a controversial subjectEI. However, there is general agreement that one important 
phase, which may actually be observed in sufficiently clean materials, is a homogeneous superconducting phase 
characterized by an energy gap vanishing only along the nodal directions (p x = ±p y in a material with tetragonal 
symmetry) and possessing conventional quasiparticle and supercurrent excitations. The low temperature behavior 
of a d-wave superconductor is described by a general action, which depends on four parameters which are defined 
below. The behavior of the parameters as the Mott phase is approached is seen to reveal information about the 
underlying physics of the Mott transition and governs the structure of vortices and thereby the doping dependence 
of such macroscopic quantities as the upper critical field. In this paper (which is largely a review) we describe the 
information which may in principle be obtained from the low T properties, summarize the present status of the 
experimental values of the different parameters, and (although the presently available data are neither fully consistent 
nor complete) outline the apparent physical implications of the results. Two subsequent papers are planned: one 
presenting the derivation of the general low temperature action from different microscopic theories and one using it 
to analyse the vortex properties (in particular iJ C 2 (T) and paraconductivity) in more detail. 

Specifically in this paper we consider a homogeneous d-wave superconducting phase, assuming in particular that 
there is no spontaneously broken time reversal symmetry and there are no other relevant excitations apart from 
the quasiparticle and phase fluctuations. (We also restrict attention to two dimensional models, but this restriction 
can easily be lifted if desired). We find that a crucial parameter is what we define as the effective charge of the 
quasiparticles. This can be determined at low T by relating the observed electronic specific heat and photoemission 
(which essentially measure the number and velocity of quasiparticles) to the temperature dependence of the London 
penetration depth (which is related to the ability of these states to carry an electrical current). Different scenarios 
for the approach to the Mott transition produce striking differences in the variation, with doping, of this charge. In 
particular, theories involving some form of spin charge separation seem to lead to a vanishing of the quasiparticle 
charge as the Mott phase is approached. We analyze available data to determine which scenario actually occurs. 

One very important feature of a superconductor is the structure of vortices introduced by thermal fluctuations or 
via a magnetic field. A vortex is characterized by a core, which may be defined in two ways: either via the density 
of quasiparticle states, which is higher near the core and drops as one moves away from the core region, or via the 
supercurrent, which varies as 1/r far from the core and drops to zero inside the core. In several of the theoretical 
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approaches to the physics of the lightly doped Mott insulator (including the one which seems to best fit the data 
discussed below), the core size as defined from the supercurrent may become very large. This, combined with the 
behavior of the quasiparticle charge, has remarkable implications for the size of the critical region, for the behavior 
of the upper critical field H C 2 and for the physics of the superconductor-insulator transition which must occur as the 
doping goes to zero. These implications were pointed out irJ3 and will be further analysed by us in a subsequent 
paper. 

The rest of this paper is organized as follows. In section II we present the low energy, long wavelength theory of a 
d-wave superconductor near a Mott insulator. In Section III we discuss the available experimental evidence concerning 
the value and doping dependence of the parameters of the theory. Section IV presents the limits of validity of the low 
energy action we discuss, along with discussion of the physics occurring when these limits are exceeded. A conclusion 
discusses the physical implications of our formalism and findings. The Appendix gives the details of calculations of 
the field dependence of the specific heat and superfluid stiffness in the vortex state. 

II. LOW ENERGY THEORY 

At low temperatures the state of all (even underdoped) superconducting cuprates seems to be a conventional d x 2_ y 2 
superconductorotm: it is described by the usual low energy degrees of freedom, namely a superconducting phase vari- 
able 4>(r, t) and, as will be discussed in more detail below, apparently conventional fermionic quasiparticle cxcitationsQ. 
Gradients of the phase correspond to supercurrents. Longitudinal supercurrents lead to charge fluctuations which are 
coupled by the long-ranged Coulomb interaction. Transverse supercurrents are described by a phase stiffness p s whose 
long-wavelength limit may be deduced from measurements of the London penetration depth. In high-T c materials, 
the low-T limit of p$ is strongly x-dependent, vanishing roughly linearly as x — > 0. This behavior is understood as 
a consequence of the suppression of charge fluctuations as the Mott insulating phase is approached, and appears to 
be related to the decrease of the resistively defined T c as x — > 0. Indeed in the very early days of high-T c Uemuta 
and co-workers showed that in underdoped materials the ratio T c (x)/p s (T — * 0,x) was essentially x-independenttl. 
At roughly the same time this behavior was shewn by a number of workers to follow naturally from theoretical mod- 
els of superconductivity near a Mott transitionO, and later Emery and Kivelson argued that the behavior could be 
understood in a more model-independent way as a consequence of the unusually small phase stiffness characteris- 
tic of high-T c materialsEl Recently Corson et. al reported direct evidence that in underdoped Bi&Sr2Cu20% the 
superconducting transition was indeed of the vortex-unbinding type driven by a small phase stiffnessEll. 

At scales less than maximal value of the gap, Ao, the physics of a two dimensional superconductor with tetragonal 
symmetry and a d x i_ y i gap function is described by an effective Lagrangian density C involving the phase (f> of the 
superconductor and quasiparticles excited out of the superconducting condensate: 

£ £(f) ~t~ fornix (1) 

Here Cp = d T — Hp is the usual Dirac action describing the 'nodal quasiparticles' excited in the vicinity of the nodes 

of the d-wave gap function. In a superconductor, the fermionic energy spectrum is given by E p — ± y 1 ' e 2 + A 2 with 

£ p the energy spectrum of the underlying fermions and A p the gap function. For a d x 2_ y 2 gap function, A p vanishes 
linearly in the four nodes, i.e. for (±7r, ±tt). It is convenient to measure momentum from the fermi point in the 
nodal direction and to parametrize the fermion dispersion by two velocities: one, vp, of the order of the underlying 
fermi velocity describing motion perpendicular to the direction in which the gap varies, and one, i>a, describing the 
opening of the gap and of order Ao /pp, obtaining 

E p = ^{v FPl ) 2 + (v A p 2 ) 2 (2) 

We take the fermions to be normal ordered in the basis which diagonalizes Hp, , so the contribution of the negative 
energy (filled Dirac sea) states is subsumed in the phase Lagrangian density which we write as 

^ = \ {di<t> + 2ieAi) * p%(r - r') * {d j( j> - 2ieA j ) (3) 

Here i is a Cartesian direction, A is the vector potential and we have allowed for non-locality in space so the * 
represents convolution. 

The quantity p 1 ^ is a diagonal matrix with dimension of energy /length 2 . Its components p x s $ (r) = p v S Q (r) = p s o(f) 
(we assume tetragonal symmetry) are related, in the absence of quasiparticle excitations, to the conventionally defined 
superfluid stiffness p s (measurable, e.g. from the London penetration depth) by 



p s (T = 0,H = 0)= / d 2 r Pso (r) 



(4) 



3 



Systems near a Mott transition are characterized by a low density of mobile charges, and we therefore expect that 
Pso( r ) has a length dependence set by this low density. 

The term C m i X gives the coupling of the phase fluctuations to the nodal quasiparticles; it may be written 

C mix = (\ d M r ) +ieA fl (r)j ■ e lq ' r Z^lf F tP+ +q/2aa iP p _ q/2taa (5) 

Here a = 1...4 labels the four nodes of the d-wave state and Z e is a phenomenological constant which we will show 
below may be thought of as the charge of a superconducting quasiparticle. It may depend on position relative to the 
node and on the proximity to the Mott transition and will be discussed in more detail below. It has been stated in the 
literature that one generically has Z e = 1; but this is now known not to be correct. Eq. C m i X is the long wavelength 
limit of a more general action involving also 'pairbreaking' terms such as r tp +r tp + with coefficients of order q. 



III. PHYSICAL CONTENT 



The low energy, long wavelength theory is described by four parameters: p$a, vp, va and the quasiparticle charge 
renormalization Z e .. To see how these parameters may be determined experimentally, we integrate out the fermions 
in the presence of static, slowly varying superflow field and vector potential, which enter via the gauge invariant 
combination 



(J = (V(j){r, t) - 2iel^ 



(6) 



and are taken to be slowly varying on the scales set by ps and the fermions. We obtain for the two dimensional free 
energy density 



F sta uc(3) = \p%Q 2 - 2T /Zj dEN ( E ) In 1 + exp[-(£ + X -Z<$ ■ t a )/T 



(7) 



where the 2 is for spin, the sum {a) is over the four nodes of the Dirac spectrum and we have introduced the single-node 
single-spin density of states per unit area 



N{E) = I -^6{E ~ E p ) = — - — 
' (2tt) 2 V PJ 2ttv f v a 



(8) 



The specific heat may be obtained by differentiating Eq. twice with respect to T and is 
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(9) 



The integral may easily be evaluated numerically for given Q, T. Analytical results exist in the limits Q/T — > and 
Q/T — * oo. The zero-field specific heat coefficient (per unit area) C(B = 0,T) is 



C(B = 0) _ 18C(3)T 

T TTVpVA 

while in the high field low-T limit we obtain (after symmetrization) 



C(Z e vQ » T) 
T 



E 



nZ e 



\2vpVA 



(10) 



(11) 



Averaging Eq [ll] over the Q-distribution characteristic of a vortex stateleads to the 'Volovik' predictiorS of a B 1 / 2 
field dependence of the specific heat if B C 2 >> B >> B c \. In principle the result depends on the nature of the 
vortex state and on the relative angle between the lattice vectors characterizing the vortex lattice (if any) and the 
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directions corresponding to the gap nodes in the superconducting state. We have evaluated the averages and find that 
the dependence is weak: 



C(B > $ v 2 F /T 2 ) _nZ e ( B 



T 3v A V$o 



1/2 

A (12) 



where A = ir/2 for a square vortex lattice with nodal direction aligned with the vortex lattice vector and A = 1.52 for 
45 degrees misalignment, similarly for a triangular vortex lattice A — 1.7 with 5% variations as the angle is varied. 
Some details of the calculation are presented in the Appendix. 

Similarly the differential penetration depth is given by differentiating F twice with respect to Q . For an arbitrary 
current distribution ps is a tensor: 

pf = Pso6 ab - £ IF^A r dx - (13) 

4ttv f v a Jo 



cosh 2 



4T 



where a, b are specific Cartesian directions_and v a is the component of vf in the a direction. Eq. ( |l6| ) describes among 
other things the nonlinear Meissner effectEj: note however the importance of the charge renormalization factor Z e . 
In the weak field limit, p s b = ps5 a b with 

\n{2)Z<*v F rr \n{2)Z^v F C{B = 0) 
PS{T) = PS0 ~ 2nv A T = PS0 ~ T— (M) 

The factor Z e is essentially the Landau parameter introduced in previous workQ. Comparison ofEqs. | and [[J 
shows why Z e is more appropriately interpreted as the quasiparticle charge renormalization. The usual f-sum-rulc 
(Ferrel-Glover-Tinkham) arguments imply that the change, with temperature, in the condensate fraction is balanced 
by an increase in the 'normal' conductivity due to quasiparticles. Now the quasiparticle conductivity is determined by 
the number of carriers (which follows from the specific heat, which essentially counts excitations) and their velocity, 
(which may be determined from photoemission) . Any remaining discrepancy with the observed dps / dT must then be 
due to their charge, i.e. to the factor Z e . . 

At large Q and low T we find the current-dependence of the superfluid stiffness to be 

P a s b (Q,T = 0) = PS0 S ab - Z* 3 « y^\^. 1?a (15) 



Wttvfva 



Calculations similar to those for the specific heat yield, for a vortex lattice with square or triangular symmetry, 
P c s = PsSab with 

PS (B,T = 0)= PSO -A-—ZJ— (16) 



IV. EXPERIMENTAL VALUES 



A. Overview 



The important parameters of the theory, ps, Z e ,vp,v\ may be determined from experiment. Of these, the crucial 
parameter is Z e . Unfortunately, the present situation is unclear because different determinations do not agree; also 
most measurements determine only combinations of the fundamental quantities, so that uncertainties in one propagate 
into uncertainties in another. In the following sub-sections section we discuss the available data for each of the three 
parameters, and then in a concluding subsection summarize the results and outstanding questions. 



B. Ps 

The T = superfluid stiffness has been reasonably well determined by muon spin rotation experimentsilll and 
decreases as the Mott insulator is approached. The decrease is apparently roughly proportional to hole doping. We 
regard this result as well established and we do not discuss it further. 
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C. v F 

Angle-resolved photoemission measurements yield t>i?0;. At present the generally accepted value fr" - i.high-T e ma- 
terials (both optimal and underdoped) along the zone diagonal and in the superconducting state isEjiij 

v F = 1.8[eV - A] (17) 

The velocity apparently increases slightly as doping is decreased. There is general agreement concerning the value 
and doping independence of the velocity (note that even undoped materials exhibit zone diagonal velocities of ap- 
proximately this magnitude). We regard this parameter as having been reasonably reliably established. 



D. V A 



The parameter va may be obtained in three ways: from photoemission, from zero-field specific heat, and from 
thermal conductivity. Each method is subject to uncertainties, as outlined below. 

Photoemission investigations of the form of the superconducting gap near th&Jiodes reveal a broadened structure, 
with a nonvanishing density of states in a small arc around the zone diagonaEJ. This could be an intrinsic effect, 
indicating a non-d-wave form of the gap function, or it could be due to pairbreaking or other sample and surface 
imperfections. However, evidence that the gapless arcs have a non-intrinsic origin is provided by penetration depth and 
thermal conductivity measurements discussed below, so we take this view here. An estimate of va from photoemission 
may be obtained by combining the gap maximum value Ao, the standard cos(20) d-wave form and the arc length 
from the zone diagonal fermi point to the gap maximum point, which is roughly n/^/2b with b the lattice constant, 
leading to 

VA = (18) 

Estimates for the gap maximum range from 30 — 40meV in optimaLl^-BCO (with the large-values in the direction 
parallel to the chains and the smaller in the direction perpendicularly to <i0meV in BSCCOcn leading to 

v A = 0.13[eV — A] (BSCCO) (19) 
v A = 0.09 - 0.12[eV - A] (YBCO) (20) 

Available photoemission evidenceH suggests that Ao and therefore va if anything increase with decreasing doping; 
suggesting (if we interpret the maximal gap observed in the (tt, 0) direction as superconducting gap) that va increases 
with decreasing doping. These estimates rely on the assumption that everywhere in the zone the observed gap has a 
superconducting origin. While this assumption has been used by many workers, and appears to be supported by the 
good argeement between the simple d-wave form and the fia^ of£3, diffferent interpretations exist in which the gap 

•^1 



in underdoped materials has a non-superconducting origii 

Eq. ( |l0| ) shows that measurements of the low temperature specific heat yield the product v f va- Because we regai 
the value of vf as reliable .-these measurements yield a value for va- In optimally doped YBCO, specific heat datat 
yield (in present notationsEJ) 

v F v A = 0.06 [eV - A} 2 (21) 

or 

va = 0.033 [eV - A] (22) 

This value is far outside the range of va suggested by photoemission. The authors of Ref@ suggest that the 
discrepancy occurs because there are additional contributions to the measured low-field specific heat (for example 
from chain states) which should not be included in the comparison between the model and data. This idea is consistent 
with recent microwave conductivity measurementscJ which find evidence for a large density of gapless excitations 
associated with the chains. An alternative possibility is that the gap function does not have the simple cos(2#) form 
often assumed, but instead is less strongly angle dependent near the nodes, so that va is not well estimated from the 
maximum gap value. Reliable measurements of the low temperature specific heat for BSCCO or underdoped YBCO 
are not available. n nn 

Thermal conductivity measurements yield values for These results rely upon a theoretical 'universal 

limit' expression for the low temperature limit of a transport coefficientEj, and upon the assumption that this low 
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temperature limit has been experimentally accessed. MeasurementsQ yield vf/va = 19 for optimally doped BSCCO. 
For YBCO a strong doping dependence is-, found. As doping is decreased the ratio drops from about 19 for a 
presumably slightly pjjerdoped FBCOe.ggaEiample to 14 for a putatively optimally doped YBCOg. 95E3 to 8 for the 
60 phase YBCOq. 73 .E3. These estimates suggest that va rapidly increases with underdoping. 

v A = 0.095 (BSCCO, overdopedYBCO) (23) 
v A = 0.13 (optimally doped YBCO) (24) 
v A = Q.2(underdopedYBCO) (25) 

These data are roughly consistent with the va inferred from the gap maximum found in photocmission; however one 
should bear in mind that the increase in v A found as doping is decreased corresponds to a decrease in the value of 
the 'universal limit' thermal conductivity. This could arise from an inhomogeneous sample (in which not all of the 
material is superconducting) or possibly from novel physics (not included in the basic action studied here) suppressing 
the ability of the quasiparticles in a doped Mott insulator to carry heat. 



E. Z e 



The crucial quantity Z e appears in combination with vf, va and so values are subject to uncertainties, particularly 
in the value of v A . 

Z 2 v 

The temperature dependence of the penetration depth yields the combination . In YBCO certainly and in 
other high-T c materials, probably, the temperature dependence of the penetration depth in the directi^it-transverse 
to the chains (if any) is only weakly material-dependent, and is linear at low T with the slope given byME3 

' 0.7 (YBC06.6,6.9) (26) 



d(k B T) 

« 0.9 (BSCCO, optimal) (27) 

Note that this linearity is inconsistent with the presence of the "gapless arcs"E3 in the electronic spectrum. From Eq. 
[T4] we then obtain 

= 6-8 (28) 

VA 

The ability to determine Z e by combining penetration depth data with values for vf and va (obtained for example 
from thermal conductivity data) was noted by Chiao, Taillefer and co-workersQ; the values obtained from the thermal 
conductivity data discussed above then yield 

Z e = 0.7 (Optimal BSCCO; OverdopedYBCO) (29) 
Z e = 0.8 (Optimally dopedYBCO) (30) 
Z e = 1 (60K YBCO) (31) 

The magnetic field dependence of the specific heat yields -^j where A is a constant (discussed above)relating to 

the current distribution in the vortex lattice. In optimally doped YBCO, high-field specific heat datan3 yield (in 
present notations^!) 

^L=0.09[eU-A] (32) 

Use of our estimate A « 1.70 and the range quoted above for va yields 

0.6 < Z e < 1 (33) 

Recent microwave conductivity measurements^ reveal an additional difficulty with the quantitative extraction of Z e 
in YBCO: in this material the deviations from tetragonal symmetry are found to lead to a strong (~ 50%) variation 
in the plane conductivity (which can be separated from the chain conductivity) between electric field parallel to the 
chain direction and antiparallel to it. This anisotropy has not been taken into account in our analysis. 
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F. Summary 



In summary, at present the experimental status of the parameter values characterizing the superconducting state 
is not completely satisfactory. The specific heat results for optimal YBCO suggest rather smaller v& values than 
are found by other measurements, and photoemission and some tunnelling data suggest that wa decreases as doping 
is reduced, while other measurements including thermal conductivity suggest that it increases. The available data 
suggest however that the crucial parameter Z e is of order unity and is only weakly dependent on doping. Particularly 
compelling in this regard is the observed weak doping dependence of dps/dT, combined with the doping independence 
of Vp, and the indications that «a increases with decreasing doping. These indications suggest that Z e is of order 
unity and if anything increases as doping is decreased. Data contradicting this conclusion exist. Further experimental 
information would be very helpful. 



V. RANGE OF APPLICABILITY OF LOW ENERGY ACTION 



A. Overview 



The results presented above constitute the leading temperature and Q dependence about the T = 0, Q — limit 
because they are nonanalytic in the standard expansion parameters (T/Eq) 2 , (vf Q/Eq) 2 where Eq is a 'microscopic' 
energy scale (for example the BCS gap amplitude Ao). We expect the expansion ceases to hold when the correction 
terms are of the order of the leading terms and in particular when the corrections to pso are of the order of pso- 

One source of correction terms are terms of the order of Q A in the phase part of the action. The usual expectation 
from study of quantum critical points is that these become important when 

Q ~ Q* = (pso/E^) 1/2 /h (34) 
~ pso (35) 

where E^,^ are 'microscopic' energy and length scales which do not vanish as the Mott phase or other critical point 
is approached. 

Another correction occurs when the fermionic terms become of the order of the leading terms, i.e. when 

T ~ T *=^fe) (36) 

or 

Roughly, if the fermionic terms determine the limits of validity of the low T, Q expansion, then the physics of the 
nonsuperconducting state is dominated by electrons and is expected to be more or less a conventional metal, whereas 
if the phase terms set the limits then fermions are irrelevant at the superconducting-non-superconducting critical 
point and the physics is presumably bosonic. 

In general, the limit of validity of the low T low B expansion signals the destruction of the superconducting state. 
We shall discuss the superconducting non-superconducting transition on the assumption that the physics is strictly 
two dimensional. While this is a reasonable approximation for high-T c materials, it is important to bear in mind 
that ultimately a crossover to three dimensional critical behavior will occur and that the parameter controlling the 
crossover is the inverse of the square of the logarithm of the superfluid stiffness anisotropy l/ln 2 (/9s||/p5_i_) which 
is not extremely small in practice, so although the two dimensional arguments provide reasonable estimates of the 
energy scales controlling T c and (as discussed below) -ff C 2, a quantitative application requires some caution. 



B. Thermal fluctuations, zero field 



In a two dimensional material the thermally driven zero-field superconducting-non-superconducting transition is a 
Kosterlitz-Thouless vortex unbinding transition. It occurs at a Tkt satisfying Tkt — ^Ps^Tkt)!^- Because thermal 
effects can only decrease ps from its T = value, the scale defined in Eq. |35| is an upper bound for this transition 
temperature. 
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Two kinds of thermal effects occur: fluctuations of the superconducting phase, and quasiparticle excitations. In the 
absence of a high density of quasiparticle excitations, longitudinal (" spin- wave" ) phase fluctuations involve unscreened 
charge fluctuations and are therefore strongly suppressed by the Coulomb interaction. In the absence of a high density 
of quasiparticles the only important excitations are vortex-antivoctex pairs, whose energetics are governed by the scale 
Pso- (In this regard we worry that the numerical studies of Ref.fEl are not quantitatively relevant to superconductors 
because these studies were based on the classical XY model, so a large contribution from 'spin-wave' longitudinal 
excitations is apparently included, whereas one would expect these to be strongly suppressed in an electronic system in 
which coulomb forces were important). Quasiparticle excitations will reduce ps{T) from its T — value. If < 1 
then the limit T* set by quasiparticle effects is more stringent. The physics of this limit (which seems to be favored 
by the data) is simply that thermal quasiparticle excitations reduce ps so that the Kosterlitz-Thouless transition 
occurs at a temperature lower by a factor of F than one would guess from pso- 



C. Field driven effects: low T 



Application of a magnetic field B > B c \ produces vortices in the superconducting order parameter. A supercon- 
ducting vortex consists of a " core" and a " far" region. In the far region superconducting excitation spectrum is only 
weakly perturbed and there is a circulating supercurrent of a magnitude jsPs/r. As the core region is approached 
the supercurrent magnitude exhibits a maximum and then decreases and the quasiparticle excitation spectrum fun- 
damentally changes. These two effects are distinct and define two core sizes, £ curr at which djs/dr = and ^ exc 
at which the excitation spectrum changes. In a conventional superconductor £ ea .--ss t; curr = up/A. In the high T c 
context, substantial attention has focused on ^ exc (which is apparently very shortu) and on the possibility that the 
change of the excitation spectrum is not simply a collapse of the superconducting gap (as in conventional materials) 
butriag tead i nvolves the appearance of a new form of long range order, for example antiferromagnetism or staggered 
fluxa'LJr 7 !! 38 !. Here we wish to focus on £ cmr which in a lightly doped Mott insulator may be much greater than £ e xc- 
Writing jg = SF/ 5Q a and using Eq leads to 
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3 S = PsoQ 
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^<(V Q -^) 2 e(V a -4) (38) 



Taking Q to be parallel to a gap node and of magnitude 1 jr we find that the current is maximal at £ curr = 
provided that £ cnrr is greater than the scale over which pso varies. Eq |3^ shows again the importance of the 
doping dependence of Z e . If (as available data seem to suggest), Z e remains constant and pso ~ x then £ cttrr ~ x~ l , 
whereas gauge-theory based modelscIJ (including, we believe, those discussed inertia) lead to Z e ~ x so that £ c urr is 
controlled by the scale dependence of ps, implying t; curr ~ x~ x / 2 . Typical numbers for optimally dopecL-i^-BCO are 
Pso ~ lOmeV, vp / Va ~ 15, implying £ cur . r [A] « I00Z 3 , roughly consistent with muon spin rotation datalLj, although 
of course uncertainties in Z lead to large uncertainties in the numerical estimates. 

The magnitude of £ curr is important because H C 2 is essentially the field at which the vortex cores overlap, and for 
the resistive transitionit is natural to use the 'current' definition of the vortex core size. Essentially this argument was 
given by Lee and WenEl who were among the first to emphasize the importance, in the high-T c context, of the scale 
over which the supercurrent varied and (on the assumption that Z e = 1) concluded that iJ c2 ~ x 2 . Future papers 
will examine in more detail the assumption that £ curr is the correct measure of the core size to use in estimating H C 2, 
but the plausibility of this claim may be seen for example from Eq. fll(i| ) which shows that when field becomes large 
enough to suppress ps by a factor of order unity the intervortex spacing is of the order of £ curr . 



VI. CONCLUSIONS 



We have presented and compared to data a general theory of low energy properties of a d-wave superconductor. 
The theory has four parameters: the T = 0, H = superfluid stiffness pso; the velocities vp and i>a describing the 
Dirac spectrum of d-wave quasiparticles, and a quantity Z e which expresses the coupling between quasiparticles and 
phase fluctuations and which we argued should be interpreted as the charge of the nodal quasiparticle. The behavior 
of these quantities contains information about the physics of the approach to the Mott transition, because different 
theoretical treatments of doped Mott insulators predict (or assume) quite different variations of these parameters with 
doping. The behavior of these quantities controls many aspects of the physics: in particular, the size (as defined from 
the supercurrent distribution) of superconducting vortices. I— ■ 

Two widely discussed theoretical approaches are the Brinkman-Rice-dynamical mean field theoryEia and the slave 
boson gauge theorycj The essential ingredient of the Brinkman-Rice theory is a self-energy £ with a strong frequency 
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dependence and a negligible momentum dependence. This leads to a Mott transition driven by a divergent dY,/doj 
implying Z e independent of x and vp ~ x. The latter prediction is in apparent contradiction to photoemission data. 
The essential assumption of the gauge theory approach is that current is carried by a small density of holes doped into 
a spin liquid environment. The fcrmionic excitations of the superconducting state are combinations of hole and spin- 
liquid states and the low density of holes leads to a small charge Z e ~ x. In other words, the quasiparticles become 
more neutral as the-,Mott insulating phase is approached. One may think of this as a precursor of the 'nodal liquid' 
phase discussed inHil. This idea also appears to be inconsistent with the available data, although, as emphasized in 
Section III the available data are not entirely consistent. Further, and perhaps most important, complete information 
is not yet available for underdoped (especially strongly underdoped) materials. We urge the experimental community 
to settle the issue of the data consistency, in order to finally establish the relevance of the Brinkman-Rice and gauge 
theory approaches to the physics of high-T c . 

Our understanding of the presently available data favors the hypothesis that Z e and vf remain constant as ps — ► 0,. 
This result would appear to rule out both the Brinkman-Rice and gauge theory descriptions of the Mott physics of 
high T c materials, and it is therefore interesting to understand the origin of the discrepancy. One common feature 
of the two approaches is that in them the Mott physics affects all of the Fermi surface in the same manner, so the 
reduction in charge stiffness is described by a reduction in velocity or in quasiparticle charge. If neither of these effects 
occurs, then the reduction in charge stiffness must be driven by a reduction in 'effective fermi surface area'. In other 
words, it seems likely that in high T c materials the crucial missing ingredient is a large, doping dependent variation 
of the parameters around the Fermi surface so that all superfluid properties arise from condensation of fcrmions in a 
narrow and doping dependent range around the nodal direction. 

Consider for example Eq. (|l^) which describes the reduction of the ps due to depairing of the nodal quasiparticles 
by a non-zero superflow. A phase gradient of magnitude Q depairs electrons in an angular range 69 ~ Z e vpQ/v^. In 
an underdoped material it seems that Z e remains of order unity while ps becomes very small. Eq. ( p"5| ) then implies 
that exciting quasiparticles in a range 69 <C tt/2 will reduce ps to zero, i.e. that all or most of the supercurrent is 
carried by electrons small patches, of angular size 89 ~ psVA/Z e v F centered on the nodal points. Within this picture 
an interesting question is the behavior of Z for angles 9 > 69. Because ps cannot become negative, the quasiparticles 
must in some manner decouple from the superfluid fluctuations (i.e Z e must become small in these regions). 

There is to our knowledge no microscopic theory of the narrow patch situation described above which is consistent 
with all data. One possibility is a commensurate long range order which opens a large, doping dependent gap around 
the antinodal points (tt, 0) which kills most of the Fermi surface leaving only hole pockets around the diagonals which 
then acquire a small amplitude superconducting gap. One example of this phenomenon would be the 'd-density 
wave' state. Another would be some form of antiferromagnetic or 'stripe' order. Two crucial consequences of such 
physics are a broken symmetry (which should be detectable in various spectroscopies) and a small va (determined by 
the observed T c ). We think that the available data do not favor this proposal. Crucially, the thermal conductivity 
measurements suggest that ua increases when doping is decreased. An alternative which is at least qualitatively 
consistent with the data is preformed (d-wave) pairs which are made mostly from the electrons near (tt, 0) regions and 
which do not contribute to any transport. For examplejfj proposed a theory in which a large mass in the (tt, 0) region 
prevented the gap maximum regions from contributing to transport. We see here that an alternative is a small Z e . 
Unfortunately there is no controlled microscopic theory which yield this physics, although uncontrolled but interesting 
extrapolations of scaling equations have been argued to lead to this physica-3. 

To summarize: elucidation of the experimental support (or lack thereof) for the 'patch picture' and (assuming it is 
relevant) clarification of its theoretical basis are two important challenges for the future. 

Acknowledgements: We thank D. Bonn, P. A. Lee, M. Chiao and L. Taillefer for very helpful discussions and L. 
Taillefer for sharing unpublished data. We acknowledge support from NSF-DMR-00081075. 



VII. 



APPENDIX: CURRENT DISTRIBUTION IN VORTEX LATTICE: H cl << H « H c2 



A. 



Formalism 



In the limit H C 2 » H » H c \ we have, for the supercurrent distribution, 



3 8 W = PsQ(r) 



(39) 



with ps the superfluid stiffness and the gradient of the superconducting phase field given by 




(40) 
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The quantity appearing in the expression for the field-induced specific heat for a vortex lattice oriented at angle 9 to 
the gap node direction is 

C{9) = -L J dxdy \Q X cos{6) + Q y sin(0)| (41) 

where the integral is over the unit cell of the vortex lattice and Ay is the area of this cell. The result has dimension 
of length" 1 . It is convenient to measure lengths in units of the inter- vortex spacing a and to normalize the result to 
the square root of the vortex density ny = B/<$>q. Thus we write 

C (9) = n v c(9) (42) 

and compute c(9) for square and triangular vortex lattices. 

B. Numerical evaluation, square vortex lattice 

We consider a square vortex lattice of lattice constant a, so the vortices sit at positions nax + may. n v = a" 2 . Eq 
gives 



Qx = £ - _?;-.T. — ( 43 ) 



1/2 (m - y/a) 

(n — x/a) 2 + (to — y/aY 

Q x {y,x) (44) 



Consider Q x . The sum over y is most conveniently evaluated in Fourier space by writing ^ m — > J dup(u) with 
p[u) = a" 1 J2 m $(v ~ ma ) = Tlk e l27rhy Substitution gives 

Qx {X,y) = ^rn]j 2 . ~ 2 , , Sm ^ 71 '^ ■ — (45) 

^ ' — ' smh (7T (n — x/a)) + sin (ny/a) 

The sum on n is rapidly convergent and may easily be evaluated numerically. 

We wish to evaluate Eq. ^l] by integrating over the region —a/2 < x,y < a/2. This is most conveniently evaluated 
numerically by inscribing a circle in the unit cell, performing the integral over the circle in polar coordinates (to 
eliminate the divergence at r — > 0) and then integrating over the remaining regions. This latter integral is over the 



region —a/2 < y < a/2; a/2 > |x| > — y 2 . We have performed this integral numerically using Mathematica; 
results are shown in the Table below. 

C. Triangular lattice 

Lattice vectors 

vi = ax (46) 
v 2 = a (^-x+^j-y^j (47) 

A general lattice vector is then nvi + tov2. The unit cell is a hexagon with area 3y/3a 2 /8. 
Eq. ^ then gives, for the component of j perpendicular to vi 

En — im — x 
w7 V (48) 

m,n (n — \m — x) +1 -^TOa — y ) 

The sum over n may again be performed-it is just the previous result with y — > x + to/2 and n — x/a — ► ^ma — y 
so that 



sinh 



'(tt (^pTO - y/aj\ + sin 2 (to7t/2 + nx/a) 
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Q x is obtained by computing the component perpendicular to a different basic lattice vector and then combining 
appropriately. 

For the Volovik effect we require J hexagon dxdy |v • j(x, y) \ . We find this is very well approximated by the integral 
over the inscribed circle. Results are shown in the Table. 



Angle 


^square 


^■triangle 





1.5708... 


1.74 


tt/8 


1.55 


1.72 


tt/4 


1.52 


1.71 



Table caption: Coefficient c defined in Eq [42] for square and triangular vortex lattice as function of angle 9 between 
lattice vector of vortex lattice and nodal direction of d-wave superconducting order parameter. 
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